On nonlinear viscoelastic deformations - a reappraisal of Fung's quasilinear viscoelastic model

This article offers a reappraisal of Fung's method for quasilinear viscoelasticity. It is shown that a number of negative features exhibited in other works, commonly attributed to the Fung approach, are merely a consequence of the way it has been applied. The approach outlined herein is shown to yield improved behaviour, and offers a straightforward scheme for solving a wide range of models. Results from the new model are contrasted with those in the literature for the case of uniaxial elongation of a bar: for an imposed stretch of an incompressible bar, and for an imposed load. In the last case, a numerical solution to a Volterra integral equation is required to obtain the results. This is achieved by a high order discretisation scheme. Finally, the stretch of a compressible viscoelastic bar is determined for two distinct materials: Horgan-Murphy and Gent

On the asymptotic properties of a canonical diffraction integral.

We introduce and study a new canonical integral, denoted I + - ε , depending on two complex parameters α 1 and α 2. It arises from the problem of wave diffraction by a quarter-plane and is heuristically constructed to capture the complex field near the tip and edges. We establish some region of analyticity of this integral in C 2 , and derive its rich asymptotic behaviour as |α 1 | and |α 2 | tend to infinity. We also study the decay properties of the function obtained from applying a specific double Cauchy integral operator to this integral. These results allow us to show that this integral shares all of the asymptotic properties expected from the key unknown function G +- arising when the quarter-plane diffraction problem is studied via a two-complex-variables Wiener-Hopf technique (see Assier & Abrahams, SIAM J. Appl. Math., in press). As a result, the integral I + - ε can be used to mimic the unknown function G +- and to build an efficient 'educated' approximation to the quarter-plane problem

One-dimensional reflection by a semi-infinite periodic row of scatterers

AbstractThree methods are described in order to solve the canonical problem of the one-dimensional reflection by a semi-infinite periodic row of identical scatterers. The exact reflection coefficient R is determined. The first method is associated with shifting the domain by a single period and subsequently considering two scatterers, one being a single scatterer and the second being the entire semi-infinite array. The second method determines the reflection coefficient RN associated with a finite array of N scatterers. The limit as N→∞ is then taken. In general RN does not converge to R in this limit, although we summarize various arguments that can be made to ensure the correct limit is achieved. The third method considers direct approaches. In particular, for point masses, the governing inhomogeneous ordinary differential equation is solved using the discrete Wiener–Hopf technique

An Efficient Semi-Analytical Scheme for Determining the Reflection of Lamb Waves in a Semi-Infinite Elastic Waveguide

The classical problem of reflection of Lamb waves from a free edge perpendicular to the centre line of an elastodynamic plate is studied. It is known that Lamb wave expansions for the displacement and stress fields poorly represent the irregular behaviour near corners, leading to the slow convergence of a series of such waves. The form of the irregularity for an elastodynamic corner is derived asymptotically, and a new solution method, which incorporates this corner behaviour analytically, is then implemented. Results are presented showing that this new approach represents the near-field and far-field behaviour very accurately, requiring very modest numbers of Lamb wave and corner modes. Further, it is revealed that the method can recover the trapped-mode phenomenon encountered in this configuration at the Lamé frequency and a specific Poisson’s ratio that we find to be approximately 0.224798

Antiplane elastic wave propagation in pre-stressed periodic structures; tuning, band gap switching and invariance

The effect of nonlinear elastic pre-stress on antiplane elastic wave propagation in a two-dimensional periodic structure is investigated. The medium consists of cylindrical annuli embedded on a periodic square lattice in a uniform host material. An identical inhomogeneous deformation is imposed in each annulus and the theory of small-on-large is used to find the incremental wave equation governing subsequent small-amplitude antiplane waves. The plane-wave-expansion method is employed in order to determine the permissable eigenfrequencies. It is found that pre-stress significantly affects the band gap structure for Mooney–Rivlin and Fung type materials, allowing stop bands to be switched on and off. However, it is also shown that for a specific class of materials, their phononic properties remain invariant under nonlinear deformation, permitting some rather interesting behaviour and leading to the possibility of phononic cloaks

Reflection from a multi-species material and its transmitted effective wavenumber.

We formally deduce closed-form expressions for the transmitted effective wavenumber of a material comprising multiple types of inclusions or particles (multi-species), dispersed in a uniform background medium. The expressions, derived here for the first time, are valid for moderate volume fractions and without restriction on the frequency. We show that the multi-species effective wavenumber is not a straightforward extension of expressions for a single species. Comparisons are drawn with state-of-the-art models in acoustics by presenting numerical results for a concrete and a water-oil emulsion in two dimensions. The limit of when one species is much smaller than the other is also discussed and we determine the background medium felt by the larger species in this limit. Surprisingly, we show that the answer is not the intuitive result predicted by self-consistent multiple scattering theories. The derivation presented here applies to the scalar wave equation with cylindrical or spherical inclusions, with any distribution of sizes, densities and wave speeds. The reflection coefficient associated with a halfspace of multi-species cylindrical inclusions is also formally derived

Characterising particulate random media from near-surface backscattering: A machine learning approach to predict particle size and concentration

To what extent can particulate random media be characterised using direct wave backscattering from a single receiver/source? Here, in a two-dimensional setting, we show using a machine learning approach that both the particle radius and concentration can be accurately measured when the boundary condition on the particles is of Dirichlet type. Although the methods we introduce could be applied to any particle type. In general backscattering is challenging to interpret for a wide range of particle concentrations, because multiple scattering cannot be ignored, except in the very dilute range. Across the concentration range from 1% to 20% we find that the mean backscattered wave field is sufficient to accurately determine the concentration of particles. However, to accurately determine the particle radius, the second moment, or average intensity, of the backscattering is necessary. We are also able to determine what is the ideal frequency range to measure a broad range of particles sizes. To get rigorous results with supervised machine learning requires a large, highly precise, dataset of backscattered waves from an infinite half-space filled with particles. We are able to create this dataset by introducing a numerical approach which accurately approximates the backscattering from an infinite half-space.EPSRC Grant EP/K033208/I and EP/R014604/

The Collision of Two Black Holes

We study the head-on collision of two equal mass, nonrotating black holes. We consider a range of cases from holes surrounded by a common horizon to holes initially separated by about $20M$, where $M$ is the mass of each hole. We determine the waveforms and energies radiated for both the $\ell = 2$ and $\ell=4$ waves resulting from the collision. In all cases studied the normal modes of the final black hole dominate the spectrum. We also estimate analytically the total gravitational radiation emitted, taking into account the tidal heating of horizons using the membrane paradigm, and other effects. For the first time we are able to compare analytic calculations, black hole perturbation theory, and strong field, nonlinear numerical calculations for this problem, and we find excellent agreement.Comment: 14 pages, 93-

A proof that multiple waves propagate in ensemble-averaged particulate materials.

Effective medium theory aims to describe a complex inhomogeneous material in terms of a few important macroscopic parameters. To characterize wave propagation through an inhomogeneous material, the most crucial parameter is the effective wavenumber. For this reason, there are many published studies on how to calculate a single effective wavenumber. Here, we present a proof that there does not exist a unique effective wavenumber; instead, there are an infinite number of such (complex) wavenumbers. We show that in most parameter regimes only a small number of these effective wavenumbers make a significant contribution to the wave field. However, to accurately calculate the reflection and transmission coefficients, a large number of the (highly attenuating) effective waves is required. For clarity, we present results for scalar (acoustic) waves for a two-dimensional material filled (over a half-space) with randomly distributed circular cylindrical inclusions. We calculate the effective medium by ensemble averaging over all possible inhomogeneities. The proof is based on the application of the Wiener-Hopf technique and makes no assumption on the wavelength, particle boundary conditions/size or volume fraction. This technique provides a simple formula for the reflection coefficient, which can be explicitly evaluated for monopole scatterers. We compare results with an alternative numerical matching method